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ABSTRACT 


This thesis investigates the feasibility of adapting the Microflown PU match vector 
sensor for underwater use. After testing the proper functioning of the sensor, the best 
materials for the capsule are detennined based on acoustic properties. The capsule is 
designed and built by NPS staff. To predict the sensitivity of the encapsulated sensor, the 
characteristic equations are modified to reflect the new medium. In order to be able to 
predict the sensitivity of the encapsulated sensor, the performance ratio is determined 
between the sensitivity in air and in the filling fluid of choice. Temperature dependency 
is introduced in the relevant parameters to be able to model the sensitivity at various 
operating temperatures. The measured sensitivity from the calibration report is then used 
to predict its performance in castor oil. The theoretical sensitivity model is verified by 
experimental data gathered from calibration studies at NUWC. The verified model is then 
used to analyze the consequences of changing critical operating parameters like the 
distance between the filaments and the operating temperature. Based on these 
calculations, recommendations are made for a better performing prototype. Eventually a 
new design is proposed that increases sensitivity significantly and is better adapted to 
operate in the filling fluid. 
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I. INTRODUCTION 

A. THE PRINCIPLE OF HOT WIRE ANEMOMETRY 

Hot wire anemometers have been used in fluid mechanics since the late 1800s. 
These thennal transducers are widely used to study turbulence and unsteady laminar 
flows. The principle of hot wire anemometry is relatively simple. A current is passed 
through a fine filament. When molecules flow by the filament they will transfer heat 
away from the wire. The rate of heat transfer will be determined by the flow speed and 
flow temperature, as well as the physical properties of the medium and filament. The 
cooling of the wire will change the resistance of the filament. The variation in the 
resistance can be measured. 

There are two basic modes of operation. The most common is the constant current 
mode, where the current though the filament is kept constant. The voltage difference to 
accomplish this will vary with the change in resistance of the filament. The voltage 
difference thus becomes a measure for the flow speed and temperature. A different mode 
of operation is the constant temperature mode. This mode will keep the filament at a 
constant temperature. In order to do this a feedback circuit will hold the resistance of the 
wire constant by varying the current through the wire to compensate for the heat transfer 
by the flow. The current thus becomes a measure for the flow speed and temperature. 

The advantage of measuring a flow by using a hot filament is the large bandwidth. 
A hot wire can measure flow speed easily from DC-flows to over 400 kHz. Furthennore a 
hot wire anemometer has very good sensitivity as well as spatial resolution. 

However, a hot wire anemometer is very easily damaged due to the exposed thin 
wires. Additionally, a hot wire anemometer is a non-linear system. This will complicate 
its calibration process [1]. Finally, a hot wire anemometer needs a minimum flow speed 
of about 1 cm/s to function properly. 
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B. THE MICROFLOWN PU MATCH 

The Dutch company Micro flown Technologies has taken a different approach in 
order to be able to measure lower particle velocities using hot filaments, using two 
parallel filaments in close proximity. The temperature upstream is slightly less than the 
temperature downstream. Therefore, the downstream wire is heated more by the upstream 
wire and vice versa when the flow changes direction. The temperature difference between 
the two wires then becomes a measure for the flow speed and direction. Using this 
principle it can operate in a flow range of 100 nm/s to 1 m/s [2], In air, this means that 
sound pressure levels (SPL) from 3.3 dB to 143 dB re 20 /./Pa can be measured. These dB 
values represent sound just over the sensitivity threshold of the human ear to well over 
the pain threshold. In seawater particle velocities are much smaller due to the increased 
impedance of the medium. 

The Micro flown Technologies “PU match” sensor adds a miniature pressure 
microphone (P) to a single axis Microflown velocity (U) sensor to fonn a package the 
size of a match. The velocity sensor measures the magnitude of the acoustic particle 
velocity along its axis. Combining the signal from the velocity sensor with a pressure 
microphone enables the sensor to resolve whether the source is in the positive or the 
negative direction along the axis and to determine the intensity of sound transmitted 
along the axis across a wide bandwidth (20 Hz-20 kHz). 

C. POTENTIAL FOR UNDERWATER SOUND DETECTION 

Although many properties are the same for both the hot wire anemometer and the 
Microflown PU match, there are two properties that motivate us only to consider the 
Microflown PU match. In water, the particle displacement and velocity are much smaller 
than in air, but the particle density is much higher. Since the Microflown PU match can 
measure much smaller flow speeds, it is much better suited for underwater sound 
detection than traditional hot wire anemometers. In addition, its beam pattern is a figure 
eight, which makes unambiguous directionality possible by adding a pressure 
microphone. 


2 



There are several reasons why it is interesting to examine the possibility to adapt 
the Microflown PU match for underwater sound detection. In the first place, it would be 
able to determine direction of sound without the need for long arrays. This is similar to 
the advantage of traditional underwater vector sensors based on accelerometer 
technology. Second, the response covers a very large bandwidth with a single sensor. 
Moreover the Microflown PU match has very good sensitivity at very low frequencies. 
Since the attenuation of low frequency sound in the ocean is much less than that of high 
frequencies, this capability could improve detection ranges. And finally, the small size of 
the sensor would enable integration in the sensor suite of very small systems. 

D. ENCAPSULATION OF THE MICROFLOWN PU MATCH 

The corrosive properties of seawater make encapsulation of the sensor necessary. 
This provides challenges common to underwater sound transducers. The encapsulation 
should be some kind of outer shell filled with a non-conducting liquid (usually oil based). 
In order to reduce reflection losses the impedance of both the shell and the liquid must be 
matched to the impedance of seawater. This impedance match will also limit refraction, 
which will introduce uncertainty in the source direction. 

It is therefore important to carefully determine the materials used for the 
encapsulation. Since resources were limited, the choice of outer shell materials was 
primarily based on the molding properties and less on acoustics. The choice of filling 
fluid was however carefully based on acoustical properties. 

This research is conducted to determine the feasibility to use this sensor for 
underwater sound detection. At this point, there is no requirement for operational 
flexibility. No effort will be put into enabling the sensor to operate in a wide depth and/or 
temperature range. 

The encapsulation approach was designed and executed by Jay Adeff (NPS staff) 
and the details will not be a part of this thesis. A general description will be provided. 
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E. STRUCTURE OF THIS THESIS 


In this thesis, the perfonnance of the sensor in air at room temperature will first be 
evaluated. This will allow us to determine if the sensor is working properly and compare 
the perfonnance underwater with its perfonnance in air, for which it was designed. The 
sensor performance underwater will then be predicted. The calibration tests to evaluate 
the performance of the sensor underwater will be executed at the Naval Undersea 
Warfare Center (NUWC) Division Newport’s Undersea Sound Reference Division and, 
due to facilities access, by NUWC staff. Data from the experiments will then be 
compared to the predictions and its perfonnance will be evaluated. The degree to which 
alterations of the cunent sensor will influence its performance underwater will then be 
examined. Finally, conclusions will be drawn and recommendations will be made for the 
way ahead. 
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II. EVALUATION OF THE MICROFLOWN PU MATCH 

IN AIR 


A. INTRODUCTION 

The Microflown PU match is a one-dimensional acoustic vector sensor. It consists 
of a conventional miniature Knowles FG-series pressure sensor and a Microflown Titan 
acoustical velocity sensor element. The sensor, as shown in Figure 1, is designed to 
operate in air between 20 Hz-20 kHz. The sensor is 45 mm long and has a width of 
3.5 mm, roughly the size of a matchstick. The upper sensor element is the conventional 
pressure microphone. The lower sensor is the Microflown Titan sensor element. In the 
sensor used for this thesis, the Titan sensor element was oriented across the stem of the 
pressure microphone element, resulting in a particle velocity axis at right angles to the 
stem. 



Figure 1. Microflown PU match (picture by Microflown Technologies) 

To verify its proper functioning and to be able to compare underwater 
performance to performance in the medium for which it was designed, the sensor was 
first evaluated in air. The primary goal was to establish the transfer function of the sensor 
and to see if the pressure and acoustical velocity sensors were coherent. The secondary 
goal was to compare the performance of the Microflown PU match sub-sensors to a well- 
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known pressure sensor; the ACO Pacific model 7046. The beam pattern of the 
Microflown Titan sensor was also detennined. These measurements were all conducted 
in the anechoic chamber at NPS. 


B. THE TRANSFER FUNCTION 

1. Method 

The method used to establish the transfer function in air was based on previous 
work by Chee Wee Ng and Jeffrey Caulk [3], [4]. To determine the transfer function, the 
sensor was placed in front of an Austin AU-15G sound source. See Appendix A for a 
comprehensive list of all equipment used in the measurements. The sound source, a guitar 
amplifier, was chosen primarily for its ability to produce loud noise over a bandwidth of 
100 Hz-2 kHz. Other acoustical properties were of less importance for this experiment. A 
calibrated ACO Pacific pressure microphone in combination with an ACO pacific 1/2” 
preamplifier was installed above the Microflown at a distance of 3.2 cm. An HP 33120A 
Arbitrary Wavefonn Generator was used to produce white noise with amplitude of 10V PP . 
The far field condition, providing the range at which wave front curvature may be 
neglected, for a frequency of 2 kHz at 20° Celsius [5] is 


FF = 


nR 1 

X 


;z-(0.032ffl) 2 

0.17m 


( 1 ) 


The distance between source and sensor was 2.04 meter and therefore well in the far 
field. 


Data acquisition during this experiment was conducted using the Mathworks 
Matlab data acquisition toolbox that enabled direct data import from the National 
Instruments cDAQ (NI cDAQ). In the system configuration used, the NI cDAQ had a 
maximum sampling rate of 51.2 kHz. The sampling rates that could be selected can be 
detennined by the following formula [6]: 

f _ f m / 256 ( ^ 

J S - -5 \ Z ) 

n 
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where fi is the sampling frequency, f m represents the master clock rate (13.1072 MHz) 
and n an integer between 1 and 31. Since the bandwidth of interest for this thesis is 
20 Hz-2.0 kHz, a chosen sampling frequency of 4267 Hz is just above the Nyquist 
frequency. This corresponds with n = 12 . 


2. The Result 


The estimate of the transfer function started by establishing the frequency range 
over which the two sensor elements were coherent. This was done using the Matlab 
command mscohere. The coherence between pressure and velocity sensor elements 
increased significantly with higher source levels. The actual transfer function estimate 
between the pressure element (p) and the velocity sensor (v) was calculated by 



(X p (k)X' v {k)) 

\x v (k)x;xk))' 


( 3 ) 


where H v (k) is the transfer function of the k lh frequency bin and X v (k) is the digital 

Fourier transform (DFT) of the time-domain signals. For better precision three datasets 
were used to compute an average transfer function. Figure 2 shows the average transfer 
function for the Microflown PU match and the coherence between the pressure element 
and the particle velocity element. The full analysis and computer code can be found in 
Appendix B. The output is coherent in the bandwidth 80 Hz-2.0 kHz. Increasing the 
source level had a significant effect on the coherence. Technical limitations of the 
equipment used prevented further exploration of the low end coherence. The transfer 
function established during the measurements is not smooth. This may be the result of 
reflections from the test setup or some airflow in the room. The amplitude of the transfer 
function from the pressure sensor element to the velocity sensor element shows the same 
trend as follows from the calibration report. The phase difference plot also follows the 
prediction from the calibration report. From this it can be concluded that the Microflown 
PU match is measuring particle vibrations and performs as stated in the calibration report. 
This will allow us to use the calibrated sensitivity as a reference value for future 
calculations. 
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AVERAGED Coherence - Hydroflown: PRESSURE to VELOCITY 
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Figure 2. Plot of transfer function for the Micro flown PU Match 


C. BEAM PATTERN 

1. Method 

To determine the beam pattern of the particle velocity element, the Micro flown 
PU match was attached to a Briiel & Kjaer Turntable system and only the output of the 
particle velocity element was recorded. A Philips compression driver was fixed to a pole 
at 2.04 m. This time the ACO Pacific reference microphone was not installed. A Stanford 
RS preamplifier was used for amplification and filtering. Pre-amplification was set to lOx 
and a 1 kHz High pass and 100 kHz Low pass filter was used to reduce noise. An Agilent 
Waveform Generator was used to produce 1000 Hz, 1500 Hz and 2000 Hz CW signals at 
1.0 V pp . 

The turntable conducted 360 degrees rotations with 6 degrees overshoot on both 
sides to ensure a continuous rotation speed during a complete circle. Data was acquired 
using National Instruments Labview and processed using Mathworks Matlab software. 


2. The Result 

During the experiment it became apparent that the vibrations of the turntable had 
a considerable effect on the experiment. Using a high pass filter at 1 kHz and measuring 
the beam pattern with frequencies above 1 kHz reduced this noise to acceptable levels. 
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The polar plot in Figure 3 shows the normalized response at 2 kHz over a 360° rotation. 
Although a frequency of 2 kHz produced the clearest results, all three frequencies showed 
the expected figure eight beam pattern provided in the data sheet. 


Beam Pattern - f=2000 Hz 
90 

1 



Figure 3. The beam pattern of the Micro flown Titan in air 
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III. SELECTION OF THE MATERIALS USED FOR 
ENCAPSULATION 

A. INTRODUCTION 

The Microflown Titan velocity sensor element uses two uncoated platinum 
filaments to measure the flow vector. Since the conducting filaments are uncoated they 
cannot be used in direct contact with seawater. To use the Microflown PU match 
underwater it therefore needs to be encapsulated. This encapsulation will affect the 
propagation of sound by reflection and refraction. Since this study is interested in both 
sound intensity and direction, it is important that the capsule has as a minimal influence 
on both these properties. 

The capsule in which the sensor will be placed consists of two basic parts; the 
outer shell and the filling fluid. For the outer shell both the material and the shape will 
determine the effect on sound transmission and direction. The filling fluid needs to be 
non-conducting and have an impedance close to that of average seawater. Other selection 
criteria were the cost and availability. Since the goal of the study is to determine the 
feasibility to develop an underwater vector sensor based on Microflown technology, it 
was considered acceptable and efficient to work with materials with less than optimal 
properties. 

The capsule was designed and produced at the NPS by J. Adeff. This chapter will 
explain the choices made in the design and some of the challenges during manufacturing. 

B. THE OUTER SHELL 

The shell was constructed using Devcon Flexane 80 urethane. This material is 
readily available and often used for hydrophones. The capsule has the shape of a cylinder 
with a hemispherical end cap (see Figure 4 and Figure 5). This guarantees a unifonn shell 
thickness to all sides. Furthermore it ensured equidistance to the center axis from all 
sides. 
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Figure 4. Capsule (designed by Jay Adeff) 
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Figure 5. Complete encapsulated sensor along cm scale 
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c. 


THE FILLING FLUID 


Castor oil, Baker DB grade, is widely used in hydrophones because of its 
excellent acoustic properties and compatibility with other transducer materials [7]. 
Furthermore, this castor oil is readily available at limited cost. Figure 6 shows how the 
characteristic impedance of castor oil (Baker DB grade) matches the characteristic 
impedance of seawater. Over the entire range of 0°C-25°C the impedance mismatch 
is <10%. 



D. SOME CRITICAL ADAPTATIONS 

The current Micro flown Titan sensor element has been designed for use in air 
only. The sensor uses platinum filaments at temperatures around 400° C to optimize 
performance. The physical properties of a liquid can be very different from a gas. For 
instance the liquid has a boiling point at which bubbles will appear. This will greatly 
decrease the sensor performance and probably reduce the life cycle. For castor oil, the 
boiling point is 313°C at atmospheric pressure. Due to increased heat transfer the Titan 
sensor element will operate at a slightly lower temperature than in air (the power 
provided to the sensor will remain the same), it will still be close to and probably over the 
boiling point of our fdling fluid. To prevent bubbles from appearing, the operating 
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temperature needs to be considerably under the boiling point. Furthermore the platinum 
used in the filaments operates as a catalyst in the dissociation of the oil. This will not only 
produce anomalies as the chemical composition of the oil changes, but it will also cause a 
honey-like residue on the filaments themselves. This can be prevented by coating the 
filaments by a thin layer of quartz [8]. 
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IV. THE PERFORMANCE UNDERWATER 


A. INTRODUCTION 

Now that composition of the prototype is known, a prediction of its perfonnance 
can be made. In order to do this some assumptions and estimates have to be made. When 
empirical data becomes available these can be validated or corrected. 

The first and most important part of predicting the perfonnance of the sensor 
underwater is to estimate how the thermodynamic properties of the filling fluid will affect 
the heat distribution. Particle displacements due to sound waves are much smaller in a 
liquid than in a gas. The density of a liquid is much greater than that of a gas, so a lot 
more molecules are available to transfer the heat. Furthermore the heat capacities of the 
two media are very different. On the other side, the diffusion speed needs to remain large 
compared to the forced convection caused by the sound wave, so that the influence of the 
particle displacement due to the sound wave can be treated as a small alteration to the 
temperature profile. The start of the prediction is therefore to calculate how all these 
factors influence the performance. 

The second part of the prediction is to estimate how the capsule will affect the 
performance. What reflection losses are expected as sound travels through the shell and 
then through the filling fluid? Is it possible to estimate the way that sound direction will 
be altered by refraction? 

This chapter will first address the influence of the filling liquid on the 
performance of the Micro flown Titan velocity sensor. Thereafter the influence of the 
packaging will be estimated. The conclusion will be an estimate of the sensor 
performance underwater. This will then be compared to empirical data in order to verify 
theory. 
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B. THE PERFORMANCE PREDICTION IN CASTOR OIL 


1. The Temperature Difference between the Two Filaments 

The sensitivity of the Microflown sensor is based on the temperature difference 
that occurs when a flow passes two heated filaments. The temperature can be calculated 
using the heat equation 

pc p (d,T + \VT)-kV 2 T = Q, (4) 

where 
p = Density 

c p = Specific heat 

T = Temperature 
k = Thennal conductivity 
Q - Generated heat 

The temperature difference will create a difference in the resistance of the wires. 
The change in resistance is therefore a measure for the flow rate. A detailed calculation of 
the temperature difference between the wires due to sound waves can be found in [9]. 
There, the term vVT (the convective term in Equation (4) ) is treated as a perturbation, 
because in air the diffusion speed is large compared to the forced convection. An 
important parameter is the heat diffusion coefficient D = k / pc p . For air at 400° C (the 

operating temperature of the Titan sensor element in air) D air « 6.5 • 10 5 nrs 1 . For castor 
oil at 313°C (the maximum operation temperature in castor oil) D « 3.2 -10 % m 2 s 1 . In the 
current sensor, with a distance between the wires (a) of 100 microns, the ratio between 
the forced convection (v) and the diffusion speed (D/a) at SPL « \40dB re 1 pPa 
compares as 
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This shows that when the sensor operates in castor oil the diffusion speed is still much 
larger than the forced convection. This implies that the convective term can still be 
treated as a perturbation if the sensor is evaluated in oil. 

Including the effect of the heat capacity of the sensors, the temperature difference 
is given by [9] 


AT = e il7tf 


2 va P 1 
D 2jzkl y if a 


x 



( 5 ) 


where: 


/ 


fo 


f = Frequency [Hz] 


fo 


D 

2na 2 


f x = 


2 nfx 

D 


D (P c r), 
271 Lh {P C P l 


Second corner frequency 


h = Thickness of the filament 


k = Thennal Conductivity Coefficient 


W 

m ■ K 


l y = Length of the wire [m] 


v = Particle velocity 


m 

s 


P = Power [W ] 


17 



a - Distance between wires [m] 

L = Width of the filament 

K x (x) = Modified Bessel function of the second kind of order 1 

The frequency / in the first tenn, as stated in the article, was replaced by the 
dimensionless frequency / , assuming that it was an error in the article. This change will 
not impact the resulting analysis of ratios as this term drops out. Calculating the ratio 
AT 

——will eliminate constant parameters from the equation, reducing uncertainty with 

^^oil 

regards to their values. Since the sensitivity is proportional to the temperature difference 
between the wires, this ratio can be used to scale the sensitivity function and thus predict 
the perfonnance. The ratio can be broken down to 
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In order to compare the performance in air to the performance in castor oil the 
calculations will be based on the same amount of acoustic intensity. The ratio between 
particle speed amplitudes ( v ) will therefore be based on the same IL. 


2. The Total Sensitivity Response 

The total frequency response of the Microflown can be approximated by [9] 


|Ar| = AT(0) 



(?) 


where 


f d = Apparent comer frequency = some arbitrary scaling factor * — , and 
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This approximation will make it much easier to model the frequency response. 
Since |A7l is proportional to the actual response of the Microflown, the performance of 

|A7l 


the sensor in castor oil can be predicted by calculating the ratio 
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Using this ratio, the expected performance of the encapsulated Microflown Titan 
sensor is shown in Figure 7. 



Frequency [Hz] 

Figure 7. The sensitivity ratio in castor oil at 313° C 

Figure 8 shows that the sensitivity ratio of the Microflown Titan element in castor 
oil at 313°C. At frequencies > 1 00 Hz. the approximation predicts less sensitivity in 
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castor oil than the analytic result of the temperature difference ratio. At 1 000Hz. the 
sensitivity ratio predicted by the approximation is twice the sensitivity predicted by the 
analytic result. 

To be able to predict the response, the calibrated sensitivity in air is needed. The 
calibration report provides a graph of the sensitivity and a formula to calculate it. In the 
preceding section the sensor response in air was confirmed to agree with this fonnula. 
Figure 8 shows the calibrated sensitivity of the Microflown PU match in air as calculated 
from this formula prior to encapsulation. 


velocity sensor sensitvity (High Gain) 



Figure 8. Sensitivity in air based on calibration report 


Since the comparison of the sensitivity in air and castor oil is calculated and the 
sensitivity in air is known, we can predict the sensitivity in castor oil by: 


Sensitivity castor oil 


Sensitivity air 
Ratio\AT\ 


(9) 


3. The Reflection Losses 

Reflection and refraction will occur due to the impedance mismatch. In the case 
of the encapsulated Microflown sensor in seawater, sound has to travel through three 
types of medium before reaching the Titan sensor element. Sound will first travel through 
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the seawater, then encounter the Flexane 80 urethane shell of the capsule. The sound then 
travels through the shell and enters the castor oil to propagate to the sensor element. At 
each boundary some portion of the sound is reflected and the direction of the sound is 
altered. Although the changes in apparent direction of sound are obviously important in a 
vector sensor, this thesis will not address this. 

The impedance of seawater is well documented and the impedance of castor oil is 
also well known. The impedance of Flexane 80 urethane is not known. Flowever, the 
shell has a very small thickness (d) compared to the wavelengths (X) on the bandwidth of 

interest. Even at 1000 Hz the ratio — < 0.01, and can therefore be considered transparent. 


The intensity transmission coefficient is given by [10] 


T, = 


4 r r 
' ' 2'1 


i r 2+ r i) 2 


( 10 ) 


where r = p x c x = characteristic acoustic impedance of the medium. 

The measurements at NUWC were conducted in seawater at 20 degrees Celsius. 
At that temperature 


' Castor oil 


1.559-10 6 ^-^, 
m 


1.446-10 6 


Pa-s 

m 


Therefore, 7) =0.999, which means that the intensity of the sound is reduced by only 

0.1% before the sound reaches the sensor element. For the comparison of the 
experimental results with the theoretical model, therefore, reflection losses can be 
neglected. 
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4 . The Predicted Performance of the Encapsulated Sensor 

Using the calibration report and the performance ratios, the predicted performance 
can be calculated. 


Predicted sensitivity, no gain 



Figure 9. The predicted performance of the encapsulated sensor 


Figure 9 displays the predicted performance of the encapsulated sensor. 
Comparing the prediction to the calibration report shows that the expected sensitivity is 
far less. Furthermore it shows that the capsule is expecting to shift the peak sensitivity to 
a lower frequency. The sensitivity of the encapsulated sensor drops off more quickly as 
the frequency increases. 

At NUWC, experiments will be conducted to determine the actual sensitivity of 
the encapsulated sensor. The results of the experiments will be compared to the predicted 
sensitivity in castor oil to verify the theory. If the theory is confirmed, it can be used to 
determine sensor adaptations to increase its performance underwater. 
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V. COMPARING THE MEASUREMENTS TO THE PREDICTED 

PERFORMANCE 


A. INTRODUCTION 

At NUWC, two sets of measurements were done using the encapsulated 
Microflown PU match. These measurements were some done a few weeks apart. During 
the first set, responses at two different orientation angles were determined, separated by 
45 degrees. First the measurements were executed along the presumed axis of maximum 
response. Thereafter, the same measurements were executed at an orientation of 45 
degrees to the assumed axis of maximum response. In the second set of measurements, 
only single orientation measurements were done. The single orientation measurements 
were intended to correspond to the maximum response axis of the particle velocity 
sensor. The aim of the measurements was to determine the responses at various Sound 
Pressure Levels (SPLs) and to determine if the sensor still exhibited its directionality 
between the two measurement angles. 

B. UNCERTAINTIES 

The nature of the sensor causes some very significant uncertainties when 
analyzing its perfonnance. First, the operating temperature is not exactly known. In air, 
the filaments are heated to approximately 400° C. The larger heat conduction in castor oil 
will reduce that temperature because the power applied will remain the same. Since the 
boiling point of castor oil is 313°C at standard atmospheric pressure, this will be 
assumed to be the temperature at which the encapsulated Microflown PU match will 
operate. Bubbles appearing at the filaments when the castor oil starts to boil will further 
influence the experimental results in an unpredictable way. 

Second, the platinum will act as a catalyst in the dissociation of the castor oil [8]. 
This process will leave a honey-like residue on the wires, reducing the sensitivity over 
time. This process is expected to influence the measurements within hours of operation. 

Finally, the exact angle of refraction is not known, because the impedance of the 

Flexane 80 shell is not known. More importantly, due to its small size, the precise 
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orientation of the Titan sensor element within the capsule is very hard to determine. Since 
the Microflown PU match has a classic figure eight beam pattern, this will influence its 
apparent sensitivity when insonified at oblique angles. 

C. THE EXPERIMENTAL SETUP 



Figure 10. Experiment setup at NUWC 


The experiments to detennine the sensitivity of the encapsulated Microflown PU 
match were conducted at NUWC by John Whitacre (NUWC staff). The most important 
part of the experiments was to detennine the sensitivity of the sensor when a source was 
located along the axis of maximum response. To determine if the encapsulated sensor 
retained its directionality, the sensor was also placed at a 45 degree angle to the initial 
orientation. During the experiments, the encapsulated sensor was submerged in an 
anechoic tank with a sound source and a well- kn own A47 SN 16 reference hydrophone. 
The bandwidth between 50-1000//zwas swept both linear and logarithmic at SPL’s 
between 90-170<ifi re 1 //Pa . Below llOdB re 1 //Pa however, the travelling field was 
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not created properly. During the experiment the temperature of the tank was 20° C. From 
the output voltage the sensitivity was calculated in dB re W / juPa . These values were 
provided for further analysis and comparison to the theory. 

The theory suggests that the current encapsulated sensor will be very insensitive 
compared to its sensitivity in air. This means that the expected output voltage is very low, 
even when high amplitude sound is directed at the sensor. The results of the first set 
suggested that the sensor was indeed very insensitive. From the first set of data it seemed 
that the sensitivity of the sensor decreased as the SPL was increased. This could indicate 
that the sensor was being overdriven. To further examine this possibility, a second, much 
more elaborate set of experiments were conducted three weeks after the first set. In this 
set the SPLs were reduced in 3 dB steps to 90 dB re 1 juPa . In total 45 measurements were 
conducted. 

D. THE RESULTS OF THE SENSITIVITY TEST 

The first test with a clear response was obtained in the first set, using a SPL of 
140dB, no gain. This was the second experiment conducted with the sensor, so the effect 
of any residue on the wires should be small. This test will also be used to determine if the 
encapsulated sensor retains its directionality. 


Sensitivity 140 dB, no gain 
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In Figure 11, the experimental sensitivity is plotted, together with the theoretical 
prediction using the approximation and using the ratio of analytic temperature 
differences. Clearly the encapsulated Microflown PU match is more sensitive to lower 
frequencies. The measurements show that the apparent sensitivity of the encapsulated 
sensor is highest around 60 Hz. It then shows an unexpected jumpy behavior in 
sensitivity throughout the measured bandwidth. On the other hand, the position of the 
peak sensitivity and the general envelope of the sensitivity seems to roughly correspond 
to our model. 

Furthermore, in this case the approximation seems to predict the sensitivity 
somewhat better than the direct calculation using the ratio of the temperature differences. 
The difference between the approximation and the exact solution is shown in [9]. That 
the approximation appears to give a better prediction of the experimental results can be 
the result of a misalignment of the sensor in the capsule. The sensor has no calibrated 
direction markings on it. The small size makes it very hard to determine the exact straight 
angle to the sensor during the encapsulation process. Because the sensor is directional, 
any misalignment will result in a reduced sensitivity in a stationary experiment. A further 
error in the assessed reflection losses caused by the impedance mismatch can further 
influence the apparent sensitivity. However, none of these factors adequately explain the 
somewhat oscillatory behavior of the sensitivity and the influence of electronic noise or 
problems with the calibration process cannot be ruled out. 
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Sensitivity 120 dB, no gain 



Figure 12. Sensitivity at 120dB, no gain 

In Figure 12, the measured sensitivity is shown using SPL of 120 dB re \/jPa 
without gain. The sensitivity over the bandwidth changes in a similar way as with the 
SPL of 140dB re 1//Pa and again the irregular sensitivity is apparent. 

Experiments conducted with SPL <H0dB re 1 /uPa gave even more irregular 
results. Sensitivity seemed to randomly vary at different SPLs. Additional experiments 
showed that the traveling field was not created well at pressure levels <110rffi re 
1 //Pa .Since the acoustic particle velocity cannot be inferred from the reference 
hydrophone when the contributions of the standing wave components are unknown, these 
experiments could not be used. 

To better understand the sensitivity of the encapsulated Microflown PU match in 
comparison with the model, all experimental data for SPL> 1 \ 0dB re 1 /uPa were 
converted to their values without any gain and then averaged. 
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0.08 


Average sensitivity, no gain 



The averaged results are shown in Figure 13, together with the predictions based 
on the approximation and the calculated temperature difference ratio. This averaging over 
eight experiments at SPLs 110-160r/fi re 1 // Pa shows a strong peak greater than the 
expected value around 60FIz followed by clear oscillatory behavior. 

This behavior could be caused by electronic noise. If 60//zand subsequent 
harmonics are the reason for the periodicity in the sensitivity, the peaks and lows would 
correspond with a harmonic frequency of 60 Hz . The initial peak is in between 56 Hz and 
64 Hz . This suggests that it is caused by electronic noise. In Table 1 the difference 
between the peaks and lows from the averaged experimental data and the 
60Hz harmonics is shown. The percentage difference seems to imply great coherence 
between the harmonics and the peaks and lows. Fiowever, the step size during these 
experiments was logarithmic. This means that the step size started with 4 Hz and 
increased over the bandwidth to 110 Hz . The peaks/lows at higher frequencies are 
therefore determined less precisely than the peaks/lows at low frequencies. To investigate 
this variation further, the raw output voltage of the measurement at 1 AOdB and 120 clB re 
\juPa were compared. Based on the 20 clB re 1 //Pa difference in SPL, a factor of 10 was 
expected in the ratio of output voltages. Fiowever, Figure 14 shows that the actual output 
voltage ratio was only about 3 and then highly fluctuating over the bandwidth. The output 
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voltage on the reference microphone showed the same fluctuations. Inspection of the 
traveling wave tube showed that, at the time of the experiments, the projector might not 
have been functioning properly. Therefore, it is justified to assume that, besides traveling 
wave components, also some standing wave components were present in the tube. The 
impact of the introduction of standing wave components is unclear since the system’s 
feedback mechanism for generating traveling waves is unknown. However, as the system 
sweeps through frequency, the large fluctuations observed on the reference hydrophone 
strongly suggest issues with the desired system response. Unfortunately time constraints 
prevented repeating the experiment in a re-calibrated tube. 
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Figure 14. Output voltage comparison 
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Further investigation to determine the cause of the periodicity in the sensitivity is needed. 


60 Hz harmonic 

peak 

low 

% difference 

60 Hz 

56.7 Hz 


6% 

120 Hz 


135.7 Hz 

12% 

180 Hz 

197.4 Hz 


9% 

300 Hz 


287 Hz 

5% 

360 Hz 

368.4 Hz 


2% 

540 Hz 


535.7 Hz 

1% 

600 Hz 

606.7 Hz 


1% 


Table 1. Comparison peaks/lows to 60FIz hannonics 

E. RESULTS OF THE DIRECTIONALITY TEST 

To determine if the encapsulated sensor retains directionality, a set of 
measurements were conducted with the sensor at a 45 degree angle relative to the initial 
angle. Although the exact refraction angle is not known, a 45 degree angle should show a 
significant drop in sensitivity at all frequencies. To determine if this was the case, all 
experiments at SPL >110r/5 re 1 juPa measured at 45 degrees were averaged and 
compared to the averaged sensitivity of the same SPLs for the initial orientation. The 
averaging was done over four data sets at both angles. 
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Figure 15. Average percentage of sensitivity loss at an angle of 45 degrees 

The results of the comparison are shown in Figure 15. Although the percentage 
loss varies over the bandwidth, it is obvious that overall the sensitivity is considerably 
lower after a rotation of 45 degrees from the presumed axis of maximum response. 
Therefore, it can be concluded that the encapsulated sensor retains directivity. The 
average loss was about 45%. The expected drop in sensitivity would be 30% from the 
measured beam pattern (see Figure 3). The difference could be due to additional 
reflection losses. An alternative cause is that the initial alignment was not precisely along 
the maximum response axis. If instead the two measurements were at angles 30 degrees 
and 75 degrees relative to the maximum response axis (and 45 degrees relative to one 
another), then this drop in sensitivity would not be unexpected. The variation of the drop 
in sensitivity, however, is significant and cannot be explained by a constant 
misalignment. This could again indicate that the traveling wave was not entirely the same 
at the position of the Micro flown during the experiments at the two orientation angles. To 
determine the actual beam pattern, additional measurements are needed, where the angle 
is changed with small increments and more data sets are gathered. 
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F. EXPERIMENTAL CONCLUSION 

The data shows that the encapsulated sensor is able to detect sound underwater. 
The envelope of the sensitivity roughly agrees with predictions using the calibrated 
sensitivity in air and the calculated sensitivity ratio or the proportional ratio in 
temperature differences between the filaments. The oscillations in the sensitivity are not 
predicted by the theory and need additional measurements to investigate. Calibration 
problems in the traveling wave tube might explain some periodicity in the apparent 
sensitivity of the Micro flown Titan sensor element. Electronic noise seems to have 
influenced measurements around 60 Hz . Further research is needed to determine its 
contribution to the apparent sensitivity at higher frequencies. The encapsulated sensor 
retains directionality. Additional measurements are needed to establish the exact beam 
pattern. 
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VI. EFFECTS OF SENSOR MODIFICATIONS 


A. INTRODUCTION 

The performance of the Microflown Titan is determined by the difference in 
temperature between the filaments as a result of a passing sound wave. A higher 
operating temperature will increase the sensitivity. As previously stated, in air the 
operating temperature of the Microflown Titan sensor element is approximately 400° C. 
The boiling point of castor oil is 313°Cat standard atmospheric pressure. Bubbles 
appearing when the liquid starts to boil will greatly reduce its perfonnance. It is therefore 
important to operate the velocity sensor element at temperatures below the boiling point. 
Flow studies using hot wire anemometry by H. Eckehnann have indicated that a 
temperature of about half the boiling point will make the effect of bubbles appearing 
insignificant [8]. Since the sensor will be used underwater, increased pressure will 
increase the boiling point and therefore also increase the optimum operating temperature. 

The higher impedance of both seawater and castor oil compared to air will result 
in a much decreased particle displacement and velocity caused by sound of the same 
Intensity Level (IL). The amount of molecules available per unit volume to carry heat 
from one filament to another however is much larger in castor oil than in air. These 
parameters translate directly into the optimal distance apart at which the filaments should 
be placed. 

In this chapter, the results of varying both the operating temperature of the 
filaments and the distance between the wires will be evaluated. This can be used to 
determine these parameters in a next prototype with the goal of improving performance. 

B. THE METHOD 

At all times the reference performance will be the Microflown Titan sensor as was 
encapsulated. This sensor, when used in air, has an operating temperature of 
approximately 400°Cand a distance between the wires of 1 OOjum . These values will not 
change throughout this chapter. 
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Several performance determining parameters are temperature dependent. The first 
step in predicting the performance of the sensor over a temperature range is therefore to 
approximate the values of these parameters over that temperature range. This was done 
by curve fitting data points using Microsoft Excel 2010 for the specific heat (c ) and 

thennal conductivity ( k ). This approximation was necessary because only sparse data 
was available. A formula to calculate the density and sound speed in castor oil was taken 
from [7]. To achieve different temperatures, the power (P) needs to become variable. 
From the power in the current sensor, the electric current (/ 0 )can be calculated. 

Assuming (/ 0 ) to be constant and using 

P = llR, (11) 

and 

R = R 0 a(AT), (12) 

we obtain 

P = lXa{AT), (13) 

with 

a =3.93-10 3 C ' = temperature coefficient, 

/„ = 9.9 -10 2 4 = current through the filament. 

The resulting formulas were integrated in Equations (5) and (7). To then calculate 
the ratios of the temperature difference and sensitivity response, Equations (6) and (8) 
were altered to be able to use a constant value for a in air and a varying value for a in 
castor oil. In Figure 16 and Figure 17 the specific heat and thennal conductivity are 
shown. Only two data points were available to determine the thermal conductivity of 
castor oil. The approximation was acceptable, however, because the total change over the 
entire temperature range was < 6% and consistent with similar oil types. 
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The approximated relationships will be treated as valid over the arbitrary 
temperature range 40°C-400°C, even though the boiling point of castor oil will make 
these approximations invalid. We can do this, because in our comparisons only the 
temperatures below the boiling point are used. 
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The temperature difference ratio now becomes 
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Using Equations (14), (15) and (9) with temperature dependency incorporated for 
the castor oil parameters, the effect of altering the operating temperature and the distance 
between the wires can be predicted. The experimental results have confirmed that the 
sensitivity in air can be scaled to the sensitivity in castor oil at a different operating 
temperature. This will be used to detennine the sensitivity of the encapsulated sensor at 
operating temperatures 40° -313°Cat any given distance between the wires. 


C. THE EFFECTS OF VARIATIONS IN OPERATING TEMPERATURE 
AND DISTANCE BETWEEN THE WIRES 

The temperature dependence of the sensitivity of the current encapsulated sensor 
is shown in Figure 18. The peak of the current sensitivity is 0.04 mV / Pa at 85Hz and at 
313°C. The sensitivity deteriorates with increasing frequency and with decreasing 
operating temperature. 
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Sensitivity prediction vs filament temperature for 100 microns 



Frequency [Hz] 

Figure 18. Temperature and frequency dependence of the sensitivity at 100 microns 

Reducing the distance between the wires will decrease the sensitivity at very low 
frequencies. Flowever, the sensitivity will be increased at higher frequencies [11]. The 
crossover point of the encapsulated sensor is < 5 Hz , so that sensitivity is in general 
increased (see Figure 19). 



Frequency [Hz] 

Figure 19. Sensitivity at 313°C and 20qm, 40 qm, 60qm, 80 qm and 100 qm 


The distance between the wires cannot be decreased indefinitely, however. The 


sensor is based on two filaments through which a current flows. As the distance between 
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the wires decreases, the induced magnetic fields and electric forces will start to obstruct 
proper performance. At what distance this will start to occur is not part of this thesis. 
However, the minimum distance will be assumed to be >20 jum . Castor oil and quartz 
(possible material for coating the platinum filaments) are very good insulators with 
dielectric constants of 4.7 and 3.8, respectively, at room temperature. 

Lowering the temperature will also lower the sensitivity over the whole frequency 
range (see Figure 18). Still, it is necessary to investigate the sensitivity at lower 
temperatures in order to prevent decreased performance by dissociation and/or boiling of 
the castor oil. The temperature at which the oil is free of boiling is dependent on the 
pressure and therefore the operating depth of the sensor. Finding the exact temperature at 
which dissociation stops is not a part of this thesis. Studies in this field by H. Eckelmann 
[8] however, suggest that below half the boiling temperature, dissociation becomes 
insignificant. For castor oil therefore, an optimum operating temperature lies somewhere 
between 100° -200°C at standard atmospheric pressure. 

The combined effect of variation of the operating temperature and the distance 
between the filaments is shown in Figures 20-22. The effects start to significantly alter 
sensitivity at < 40 jum between the filaments. 


Sensitivity prediction vs filament temperature for 40 microns 



Frequency [Hz] 

Figure 20. Temperature and frequency dependence of the sensitivity at 40 microns 
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Sensitivity prediction vs filament temperature for 30 microns 



Frequency [Hz] 

Figure 21. Temperature and frequency dependence of the sensitivity at 30 microns 


Sensitivity prediction vs filament temperature for 20 microns 



Frequency [Hz] 

Figure 22. Temperature and frequency dependence of the sensitivity at 20 microns 

For 313°C and 200° C the effect of changing the distance between the wires is shown in 
Figures 23 and 24, respectively. 
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Predicted sensitivity, no gain at 313 deg C 



Figure 23. Distance dependence of the sensitivity at 313°C 


Predicted sensitivity, no gain at 200deg C 



Frequency [Hz] 


Figure 24. Distance dependence of the sensitivity at 200° C 


From Figures 23 and 24 it can be concluded that decreasing the distance between 
the wires has a much larger effect than changing the temperature in the ranges of interest. 
Based on ratio of the temperature differences between the wires, decreasing the distance 
between the filaments to 40 jum will result in a 8dB gain in sensitivity at both 313°C and 
200°C. Decreasing the distance to 20//mwill result in a 14r/figain for both 
temperatures. A filling fluid with similar impedance and with a much higher boiling point 
might improve this gain even further. 
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According to the theory, decreasing the distance between the filaments will increase 
sensitivity rapidly as the distance between the wires gets very small. This raises the 
question as to whether the theory still holds at distances <30 jum or whether other 
influences become dominant. If the theory does hold, it would be useful to establish the 
minimum distance apart that can be achieved. 

The sensitivity peaks just below 1 00Hz, in the region where the spikes were 
observed during the experiments with the current sensor. If the sensor is already sensitive 
in this region, reducing the distance between the wires might increase sensitivity to the 
point that the sensor is overdriven. In that case, a design with two sensor elements could 
be considered. One sensor element optimized for frequencies < 150 Hz and one optimized 
for higher frequencies. Since the sensitivity peak does not shift very far as the distance 
between the wires decreases it is not necessary to add more sensor elements. The size of 
the elements (only a few millimeters) enables stacking them without any significant 
increase in overall sensor size. 

D. LIMITS IN FILAMENT PROXIMITY 

Although theory suggests that a reduction of the distance between the wires will 
increase the sensitivity significantly, technological restrictions will limit the extent to 
which the wires can be brought together. J.W. van Honschoten has experimented in air 
with filaments set30///n apart [9]. It therefore seems reasonable to assume that sensors 
with a distance between the filaments of about 30 jum can be produced. 

E. SENSOR MODIFICATION CONCLUSIONS 

Reducing the distance between the wires increases the sensitivity significantly. 
Reducing the operating temperature will have a relatively small negative effect on the 
sensitivity. Reducing the operating temperature in castor oil and coating the filaments is 
necessary to avoid boiling and dissociation of the castor oil. Literature indicates that this 
effect is negligible when the operating temperature is reduced to about 200° C. The effect 
of reducing the distance between the wires increases as the filaments are closer together. 
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This encourages experiments to determine the validity of the underlying theory at these 
distances. In future designs of the sensor a stacking of two elements could be considered 
if the sensor is found to be too sensitive at low frequencies. 
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VII. CONCLUSION AND RECOMMENDATIONS 


A. CONCLUSIONS 

It seems feasible to adapt the Microflown PU match for underwater use. 
Experiments seem to indicate that the theories to predict its performance hold when the 
sensor is submerged in oil. However, the oscillations of the measured sensitivity and the 
high peak in measured sensitivity around 60 Hz are not predicted by our model and need 
further investigation. The influence of the calibration error in the experimental set up and 
the influence of electronic noise need to be determined. 

Based on its acoustic and electric properties, castor oil is a good choice as a fdling 
fluid. Additional gain can be achieved by finding a more stable filling fluid with similar 
or better acoustic and electric properties, but a much higher boiling point. Operating the 
sensor in a liquid instead of a gas will require some adaptations to the sensor. The 
operating temperature needs to be well below the boiling point of the oil. Coating the 
filaments will be necessary to overcome catalytic effects of the platinum wires. Further 
strengthening of the wires may be required to prevent damage when operating in a much 
denser medium. 

At the current distance between the filaments, sensitivity underwater is low. 
Decreasing the distance between the filaments will increase the sensitivity over the whole 
bandwidth of interest. A gain of 14-21 dB could be achieved if the distance between the 
wires is reduced to 20-40 jum. Reducing the operating temperature to 200° C does not 
significantly reduce sensitivity at these distances. 

The encapsulated sensor retains a level of directivity. The exact beam pattern 
needs to be detennined in future research. 

B. RECOMMENDATIONS FOR FUTURE DESIGN 

The next prototype should have coated filaments to prevent the platinum from 
acting as a catalyst in the dissociation of the oil. The operating temperature needs to be 
well below the boiling point of the inert liquid of choice. Based on its acoustic properties 
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castor oil is a good filling fluid for the capsule. In this case an operating temperature of 
200° C is advised. The wire distance should be reduced to <40 jum in order to improve 
sensitivity. 


C. RECOMMENDATIONS FOR FUTURE RESEARCH 

Future research is necessary to investigate the oscillations in the sensitivity. The 
conducted experiments should be repeated in a re-calibrated traveling wave tube to 
determine if standing wave components are the cause of the periodicity in the sensitivity. 
Additional experiments are necessary to determine the influence of electronic noise. 
Subsequent research is necessary to determine if the beam pattern of the encapsulated 
sensor maintains its proper response. Finally, it needs to be determined if the current 
prediction model will hold at very small distances between the wires. 
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APPENDIX A. LIST OF EQUIPMENT 


Equipment and settings used during all measurements 

Microflown PU match consisting of: 

- Microflown Titan acoustical particle 

velocity sensor 

- Knowles FG series miniature sound 

pressure transducer 

Frequency range: 

- Velocity sensor: 0.1 Hz-20 kHz 

- Pressure sensor: 20 Hz-20 kHz 

Microflown Technologies MFSC-2 

2 channel signal conditioner 

- Gain set to ‘High’ 

- Correction mode set to ‘Off 

Post-processing: 

Mathworks Matlab Version 2010a 


Additional equipment and settings used to determine beam pattern 

Agilent type 33220A 

Arbitrary Waveform Generator 

- 1000 Hz, 1500 Hz, 2000 Hz CW 

- Vppsetto ‘1.0 V’ 

Briiel & Kjaer Turntable system type 9640 


Stanford RS preamplifier model SR560 

Pre-amplification: lOx 

High pass filter: 1 kHz 

Low pass filter: 100 kHz 

Philips compression driver 


Data acquisition software: 

National Instruments Labview 2011 

Sampling frequency: 

25 Hz 

Additional equipment and software used to detennine transfer function 

National Instruments Compact DAQ USB 

chassis model NI 9172 


National Instruments Sound and Vibration 

DAQ module model NI 9234 (2 total) 
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Austin AU-15G Amplifier 

- Overdrive set to ‘ON’ 

- Volume set to ‘MAX’ 

Hewlett Packard Function Generator model 

33120A 

- Mode: White Noise 

-Amplitude: 10V PP 

AGO Pacific Calibrated Pressure 

microphone model 7046 


AGO Pacific V” pre-amplifier model 4012 


Data acquisition software: 

- Mathworks Matlab 2011 

- Data acquisition toolbox 


Additional equipment and software analyze data NUWC 

Microsoft Excel 2010 
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APPENDIX B. EVALUATION OF THE MICROFLOWN IN AIR 


A. THE FULL EVALUATION 

The proper functioning of the Microflown PU Match sensor and all its sensor 
elements was confirmed. To be able to do this, first the transfer function was established 
between the pressure element and the velocity element. Also the transfer function with an 
independent and well kn own ACO Pacific model 7046 microphone was established and 
compared to the initial results. Finally the Noise Spectrum Level (NSL) was determined. 
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THE MATLAB CODE FOR THE TRANSFER FUNCTION 


%% Hydroflown Transfer Function 

% Program to read in Hydroflown data and compute transfer function between 
% pressure sensor and velocity sensors 


%% Clean up 
clear all 
cl c 

close all 

%% Setting up F FT 
maxchan = 6; 
numpass = 250; 

N = 1024; 

M = numpass * N; 

Fs = 4267; 

NF FT = 2 A 10; 

over lap = NFFT/ 2; 

wi ndow = hanni ng( NFFT) ; 


%maximum number of data channels captured 
%number of loop iterations in vi sampling routine 
%number of samples per channel per pass 
%t o t a I sample length per channel 
%s a mpI i n g f r e q 


%% Channel 
p r s 1 = 1; 
bl ul = 2; 
gr nl = 3; 
pr s R = 4; 
p r s 2 = 5; 
bl u2 = 6; 


Desi gnation 
% sensor 324 pressure 
% sensor 324 blu veloci 
% unused 
% ACO pressure 
% unused 
% unused 


ty 


Srefl = prsl; % hydroflown Pressure Channel 
S r e f 2 = bl ul; % hydroflown Velocity Channel 


pi ot sensor =' Mi crofl own' ; 
pi ot ref 1=' PRESSURE’ ; 
pi ot ref 2=' VELOCI TY’ ; 
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%% Loadi ng Datasets 


bindata = zeros( M, maxchan); ^initialize matrix to allocate memory 


datasets = 3; %number of data sets for averaging 


for ds =1: dat as et s %c y c I e through number of data sets 


if ds == 1 
bindata = 

i mport data( ' C: \ Users \ Marni x \ Documents \ MAT L A B \ Thesi s \ measurement s \ Mar20_ WN _10Vpp_max 
_l_60s. mat’ ); 

else if ds == 2 

bindata = 

i mportdataf' C:\Users\Marni x\Documents \ MAT LAB\Thesi s\ me as ur ement s\ Mar 20_ WN_ 10 V p p _ max 
_2_60s. mat’ ); 


else if ds == 3 


bindata = 

i mport data{' C: \ Users \ Marni x \ Documents \ MAT L A B \ Thesi s\measurements\Mar20_ WN _10Vpp_max 
_ 3_ 6 0s. mat 1 ); 


end 


end 


end 


bindata = bindata’; %change to row data 


%% C o mp u t a t i o n for Transfer Function 
channel 

g a mma _ p 2 v ( : , d s) = 

mscoherefbi ndatafSrefl, ; ), bi ndata(Sref2, 

[ P2V(; , ds), F] = cpsd(bi ndata(Sref 1,: 

[ V(: , ds), F] = pweI c h( bi ndat a(Sref2, : 
[P(:,ds),F] = pwe I c h ( bi nd a t a (Sr ef 1, : 
Hp 2 v ( ; , ds) = P 2 V ( : , ds) . / V( : , ds); ”/ 


for Sound Probe Pressure to velocity 

), wi ndow, over I ap, NFFT, Fs); 

), bi ndata(Sref2, ; ), wi ndow, overl ap, NFFT, Fs); 
), wi ndow, over I ap, NFFT, Fs) ; 

), wi ndow, over I ap, NFFT, Fs) ; 

PRESSURE to Vel oci ty 
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%% Computation for ACO Reference to Pressure element channel 
gamma_pACO( : , ds) = 

ms coher e( bi n d a t a ( prsR, : ), bi n d a t a (Sref 1, : ), wi ndow, overl ap, NFFT, Fs); 
[ P2AC0( : , ds), F] = 

cpsdjbi ndataf prsR, : ), bi ndatafSrefl, : ), wi ndow, overl ap, NFFT, Fs); 

[P(:,ds),F] = pwel ch( bi ndataf Sref 1, ; ), wi ndow, overl ap, NFFT, Fs); 

F2x(;,ds) = P2AC0(;,ds)./P(:,ds); % ACO t o PRESSURE 


end 


%% PLOT Transfer functions 

% VELOCITY element to PRESSURE element Transfer Function 
f i gu r e( 1) 
subpl ot (3, 1, 1) 
hoi d on 

pi ot (F, gamma_p2v( ; , 1), ' b' ) 
yl abel (' Coherence' ) 

t i 11 e( [' Coher ence • .plotsensor, .plotrefl, to 
, pi o t r e f 2 j, ' Font Wei g h t',' bold’) 

axi s ([ 0 2 0 0 0 0 1.1 ]) 

subpl ot (3, 1, 2) 
hoi d on 

pi ot(F, abs( Hp2v( : , 1)), ' b 1 ) 

yl abel (' Tr ansf er Function Amplitude') 

titleU'Amplitude of Transfer Function ■ '.plotsensor, .plotrefl, 
‘ , p I o t r e f 2 ], 'Font We ight', 'bold') 

axi s([ 0 2 0 0 0 0 2 ]) 

subpl ot (3, 1, 3) 
hoi d on 

pi o t (F, rad2deg( angl e (H p 2 v ( ; , 1))), ' b' ) 

t i 11 e ([ ' P h a s e Difference - .plotsensor, .plotrefl,' to 
, pi o t r e f 2 j,' Font We ight', 'bold') 

yl abel (' Phase Difference (deg)’) 


t o 
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xl abel( ' Frequency (Hz)’ ) 

axi s([ 0 2 0 0 0 - 1 8 0 1 8 0 ]) 

% PLOT averaged VELOCITY element to PRESSURE element Function 
fi gure(2) 
s ubpl ot (3,1, 1) 
hoi d on 

gamma_p2v=sum( gamma_p2v, 2) / 3; 

pi ot (F, gamma_p2v, ' b' ) 

yl abel (' Coherence' , ' FontSi ze' , 14) 

t i 11 e ([' AVERAGE D Coherence - , pi otsensor, .plotrefl,' to 

, pi o t r e f 2 ], ' FontSi z e' , 18,' Font Wei g h t ’ , ' bol d ’) 

axi s ([ 0 2 0 0 0 0 1.1 ]) 

subpl ot (3, 1, 2) 
hoi d on 

Hp2v = sum] Hp2v, 2)/3; 
plot(F,abs(Hp2v),'b') 
yl abel (' A mp I i t u d e' , ' FontSi z e' , 14) 

t i 11 e ([' AVE RAGE D Amplitude of Transfer Function - ' , pi otsensor, ,plotrefl, 
to ' , pi otref 2], ’ FontSi ze' , 18, ' Font We ight'.'bold’) 

axi s([ 0 2 0 0 0 0 2 ]) 

subpl ot (3, 1, 3) 

pi o t (F, rad2deg( angl e ( H p 2 v)), ' b’ ) 

t i 11 e ([' AVE RAGE D Phase Difference - .plotsensor, .plotrefl,' to 
', pi o t r e f 2 ], 1 FontSi ze’, 18, 'Font We ight'.'bold') 

yl abel (' Phase Difference (deg)' , ' Font Si z e ’, 14) 

xl abel ( ' Frequency (Hz) ’ , ’ FontSi ze’, 14' ) 

axi s([ 0 2 0 0 0 - 1 8 0 1 8 0 )) 

% PLOT PRESSURE element to ACO microphone Transfer Function 
f i gur e(3) 
subpl ot (3, 1, 1) 
hoi d on 
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pi ot (F, gamma_pACO( : , 1), ' b' ) 
yl abel{' Coherence' ) 

t i 11 e( [' Coher enc e - , pi otsensor, ' : ACO Microphone to ‘ 

plotrefl], 'Font We ight’, 'bold') 

axi s([ 0 2 0 0 0 0 1.1 ]) 

subpl ot (3, 1, 2) 
hoi d on 

pi ot{ F, abs( F2x(:, 1)), ' b* ) 

yl abel (' Transfer Function Amplitude') 

t i 11 e{ [' Ampl i t ude of Transfer Function - ' , pi otsensor, : ACO Microphone to 
plotrefl], 'Font We ight’, 'bold') 

axi s([ 0 2 0 0 0 0 2 ]) 

subpl ot ( 3, 1, 3) 
hoi d on 

pi o t (F, r ad2deg( angl e (F 2 x ( : , 1))), ' b’ ) 

t i 11 e( [ 1 Phase Difference - ,plotsensor, ACO Microphone to ' 

plotrefl], Font We ight 1 , 'bold') 

yl abel (' Phase Di f f er ence (deg)') 

xl abel (' Frequency (Hz)' ) 

axi s([ 0 2 0 0 0 ■ 1 8 0 1 8 0 ]) 

% PLOT averaged PRESSURE element to ACO Transfer Function 
f i gu r e(4) 
subpl ot (3, 1, 1) 
hoi d on 

gamma_pACO=sum( gamma_pACO, 2) / 3; 
pi ot (F, gamma_pACO, ' b’ ) 
yl abel (' Coherence' ) 

t i 11 e( [' AVERAGED Coherence - ,plotsensor, ACO Microphone to 1 

plotrefl], 'Font We ight', 'bold') 

axi s([ 0 2 0 0 0 0 1.1 ]) 

subpl ot (3, 1, 2) 
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hoi d on 


F 2 x = s u m( F 2 x, 2) / 3; 
plot(F,abs(F2x),'b ) 
ylabelf'Transfer Function Amplitude’) 

t i 11 e ([ ' AVERAGE D Amplitude of Transfer Function ■ pi ot sensor, ACO 
Microphone to ' , pi otrefl],' Font Weight’ bold') 

axi s([ 0 2 0 0 0 0 2 ]) 

subpl ot (3, 1, 3) 
hoi d on 

pi ot (F, r ad2deg(angl e( F2x)), ' b' ) 

t i 11 e ([' AVE RAGE D Phase Difference - .plotsensor, ACO Microphone to ' 

, pi otrefl], ' Font We i ght' , ' bol d’ ) 

yl abel (' Phase Di f f er ence (deg)') 

xl abel ( 'Frequency (Hz)’ ) 

axi s([ 0 2 0 0 0 - 1 8 0 1 8 0 ] ) 

% PLOT VELOCITY element to ACO microphone transfer function 
f i gu r e(5) 
subpl ot (2, 1, 1) 

F H = H p 2 v. * F 2x; %t r a n s f e r function to correct velocity to ACO pressure direct 

plot(F,abs(FH),’b') 

yl abel {'Transfer Function Amplitude') 

t i 11 e( [' Ampl i t ude of Transfer Function • ', pi ot sens or, ’ : ACO Microphone to 
', pi o t r e f 2 ], ' Font Wei ght'bold') 

axi s([ 0 2 0 0 0 0 3 ]) 

s ubpl ot (2, 1,2) 

pi ot (F, r a d 2 d e g (angle(FH)),'b') 
yl abel (' Phase Di f f er ence (deg)' ) 
xl abel( ' Frequency (Hz)’ ) 

t i 11 e( [' Phase difference of Transfer Function ■ , pi otsensor, ' : ACO Microphone 

to ' , pi o t r e f 2 ], ’ F o n t We i ght' , ' bol d ’ ) 

axi s([ 0 2 0 0 0 - 1 8 0 1 8 0 ]) 

%% Determine Noise Spectrum Level 
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bindata = zeros( M, maxchan); %i n i t i a I i z e matrix to allocate memory 

bi ndat a = 

i mport data( ' C: \ Users \ Marni x \ Documents \ MAT L A B \ Thesi s\measurements\Mar20_ WN _10Vpp_max 
_3_60s, mat' ); %i n i t i a I i z e ma t r i x to allocate me mo r y 

m=0. 0 5 4 3 3; %aco sensitivity 

bindata = bindata'; %change to row data 

noi se_si gnal = bindata(Sref2,:); 

noi se_si gnal = noi se_si gnal - mean( noi se_si gnal ); 

pent = p we I ch( noi s e _ s i gnal , wi ndow, overl ap, NFFT, Fs); 

pe nf =pe nf. * (F H. * c o n j (F H) ) / m A 2* 1E 1 2/ 4 0 0; %pressure equivalent noise floor 
computation re 20uPa 

pent = 10*Iog 10(penf); 

ACO = bi ndataf prsR, : ) ; 

ACO = ACO • mean) ACO); 

[PACO, F] = pwel ch(ACO, wi ndow, overl ap, NFFT, Fs) ; 

PACO = PACO/m A 2; % ACO corrected to Pa 
PACO = PACO* 1E 1 2/ 400; % ACO signal in 20uPa 

PACO = 10*1og10( PACO); 
di f f =PACO- pent; 


%% Plot Noise Sped r urn Level s 

% PLOT Noise Spectrum Level measured by ACO microphone and PRESSURE element 
f i gur e(6) 
s ubpl ot (3, 1, 1) 
pi ot(F, pent, ' b' ) 
xl abel ( ' Frequency (Hz)' ) 
yl abel( ' NSL re 20uPa A 2/Hz’ ) 

t i 11 e( [' Pressure Equivalent Noise Spectrum Level of .plotsensor, 
plotref2], ’Font We ight’,' bold’) 

axis([0 2 0 0 0 0 5 5 ]) 
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s ubpl ot (3, 1,2) 
pi o t ( F, PACO, ' r' ) 
xl abel (' Frequency ( Hz)' ) 
yl abel ( ' NSL re 20uPa A 2/ Hz' ) 

ti tl e(‘ Pressure Equivalent Noise Sped r um Level of ACO 
Mi crophone' , 1 Font Wei g h t ’ , ’ bol cT ) 

axis([0 2 0 0 0 0 5 5 ]) 


s ubpl ot (3, 1, 3) 
pi ot(F, di ff, b' ) 
hoi d on 

% p I o t (F, p e n f, ' r' ) 
xl abel ( 1 Frequency (Hz)’) 
yl abel ('NSL re 20uPa A 2/ Hz' ) 

title(' Difference in Pressure Equivalent NSL of ACO Microphone and 
Velocity','Font We ight'/bold') 

axis([ 0 2 0 0 0 - 1 1 0 ]) 

MATLAB CODE FOR DETERMINING THE BEAM PATTERN 

%% Cal c u 11 e beampattern on Microflown particle velocity sensor 
% frequency 1500 Hz 
% S a mp I i n g period 4. 0 0 e - 2 s 
cl ear; 
cl c 


c 1 c =xI s r e a d ( 2000_cropped_cl ean, xl 

tt =4. OOe-2*1engt h( clc) ; 
c1c(: , 1) =clc( : , 1) - c1c(1, 1); 
a mp =c 1 c ( ; , 2); 
dbamp=20. *1 ogl0( amp) ; 
dbamp=dbamp- mi n( dbamp) ; 
dbamp=dbamp / max] dbamp); 
angamp] ; , 1) =c 1 c(; , 1) *((2* pi ) /11); 
angamp] ; , 2) =d b a mp; 


sx' ); 

% t o t a I ti me [ s ] 

% st art at t =0 s 
% define a mp [ V ] 

% convert to decibels 
% s e t mi n i mu m t o 0 d B 

% convert time to radians 
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pol ar( angamp( : , 1), angamp( : , 2)); % plot beampattern 

titlef'Beam Pattern ■ f =2 0 0 0 Hz'); 
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APPENDIX C. MATLAB CODE TO DETERMINE ACOUSTICAL 
PROPERTIES OF CASTOR OIL 


%% Rho_0, c Calculation of Castor oil 

% Intended to gain insight in the acoustic properties of castor oil and seawater as 
a 

% function of temperture 
%% Clean up 
clear 
cl c 

close all 
%% p a r a me t e r s 

t =0: 0. 1:2 5; %Temperture range of the oil indegCelcius 

p =2 0.265; % Pressure in dBar 

SA=35*ones( 1, I engt h( t)); % Salinity [g/kg] (constant in this case) 

%% Calculation density and speed of sound Castor oil 

r ho=1000. / ( 1. 02714+7. 04e - 4. *t +9. 66e - 7. *t . ''2+3. Oe- 9. *t. "3-4. 91e - 4*( p/ 100) - 2. 633e- 
6*(p/100).*t-4.04e-9*(p/100).*t. A 2-8.8e-ll*(p/100).*t.''3+1.471e-6*(p/100) A 2+9.2e- 
9*( p/ 100) A 2. *t. A 2-3. 6 3 e - 9 * ( p/100) A 3-l. 6 6 e - 11* ( p/ 100) A 3. *t); 

c =1 5 7 0* ( 1. 0 0 0- 2. 1 5 e - 3. *t +4. Oe - 6. *t. A 2+2. 5e - 6*( p/ 100) . *t +2. 2 2 e - 3*( p/ 100) - 3. Oe- 

6* ( p/ 10 0) A 2) ; 

%% Calculation density and speed of sound sea water 
[ r ho_s] = gsw_r ho_t_exact (SA, t, p) ; 

[c _ s]= gsw_sound_speed_t_exact (SA, t, p); 

%% r h o * c products 

rho_c_c=rho.*c; % Characteristic! mpedance 

Castor oil 

rho_c_s=rho_s. *c_s; % Char act er i st i c Impedance 

Seawater 

del ta=( (abs( rho_c_c - rho_c_s)). / rho_c_s)*100; % Difference in Impedance [ %] 

%% Plot 

figure (1) 

subplot (2,1,1) 

pi ot (t, r ho) 

hoi d on 

grid on 
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pi ot(t, r ho_s, ' r ’ ) 

xl a bel (' Temper t ur e [degrees Cel c i us]' ) 
yl abel (' Densi ty [kg/ m A 3]' ) 

hi e g 1 = I egendf 1 castor oi I ' , ’ seawater' , ' Locati on' , ' South West' ); 

t i 11 e( ' Dens i t y and sound speed in Castor Oil and Seawater at 
10m' ,' FontWei ght’ , ’ bol d' ) 

subplot (2,1,2) 

pi ot(t, c) 

hoi d on 

grid on 

pi ot(t, c_s, 1 r' ) 

xl a bel (' Temper t ur e [degrees Cel c i us] ' ) 
yl abel (' speed of sound [m/s]') 

hi egl=l egend( ' castor oil','sea wa ter’,’Location’, ’South’); 

f i gur e(2) 

subplot (2,1,1) 

plot(t,rho_c_c) 

hoi d on 

grid on 

pi ot(t, r ho_c_s, ' r ’ ) 

xlabel('Temperture [degrees Cel ci us]' , ’ FontSi ze’ , 16, ' FontWei ght’ , ' demi ' ) 
yl abel (‘ I mpedance [ Pa*s/m] ’ , ’ FontSi ze’ , 16, ' FontWei ght' , ' demi ‘ ) 
hieg1=I egendf 'castor 

oi I ' , ' seawater' , ' Locati on' , ’ South We st’ , ' FontSi ze’ , 16, ' Font We ight'.’demi'); 

t i 11 e(' Char act er I st i c Impedances of Castor Oil and Seawater at 
10 m' ,' FontSi ze' , 20, ' Font Wei ght' , ' bol d’ ) 

subpl ot (2, 1,2) 

pi ot (t, delta) 

grid on 

xlabel(‘Temperture [degrees Cel ci us]' , ' FontSi ze' , 16, ’ FontWei ght' ,' demi ’ ) 
yl abel (' Difference in Impedance [%]’ , ' FontSi ze’, 16, ' FontWei ght’ , ’ demi ' ) 
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APPENDIX D. MATLAB CODE FOR PERFORMANCE 
PREDICTION 

%% Sensitivity prediction 

% Script determines the temperature differences between the wires for both air and 
Castor oil based on the 

% the formulae in the article by J.W. Honsbroek, "Analytic model of a two-wire 
t h e r maI 

% sensor for flow and sound measurements).' 1 It also mo dels the 
% approximation given in the articles. All parameters are made temperature 
% dependent to predict performance over a wide range of temperatures. 

%% Clean up 
cl ear; 
cl c; 

close all 

f o r ma t long 

%% Basic p a r a me t e r s 

f =0: 1; 1 0 0 0; 

T =313; 1; 6 7 3; 

TT =4 0: 1:4 0 0; 

TTT =6 7 3; 

p =1 0. 1 3 2 5; 

st ep=20e- 6; 

v_a =0. 0 2 4 2 0; 
air) [ ml s ] 

v_c =6. 186e- 06; 

[ m/ s ] 

[ X Y]=mes hgri d(f, T); 

% Mi c r of I own 
I y =0. 001; 

L =2e- 6; 

A =2 e - 10; 

I 0=9. 8 71 e - 2; 
a_a=100e-6; 


% frequency range [Hz] 

% T e mp e r a t u r e range [K] 

% T e mp e r a t u r e range in Celcius 

% a p p r o x i ma t e operating temperature (air) [deg C] 

% At mo spheric Pressure in d B a r 
% step distance between wire [ m] 

% a b i t r a r y particle velocity (IL 78 dB re le-12 W/ m A 2 in 

% Particle velocity (IL 78 dB re le-12 W/ m A 2 in Castor Oil) 


% Length of the f i I a me n t s [ m] 

% E s t i ma t e d half width of the f i I a me n t [ m] 

% Crosssection of a filament 
% Current based on P =12mW i n air [A] 

% distance between wire in the microflown Titan [ m] 
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a_c=20e- 6; % initial (variable) distance between wire potential 

hydr of I own [ m] 

a_max=100e - 6; % maxi mum di st ance between wires 

ns=a_max / st ep; % Number of steps 

h_pt =100e- 9; % thickness of pi at I num I ayer [ m] 

al pha=3. 93e - 3; % Temperature coefficient of resistivity Platinum [C A -1] 

h_Si Ni =200e-9; % thickness of SiNi layer [m] 

r hoc_pt =2. S5e6; % density times specific heat Pt [J/(m A 3*K)] 

rhoc_Si Ni =1.66e6; % density times specific heat SiNi [J/(m A 3*K)] 

%% P r e - a I I o c a t e ma t r i c e s 

rati o_FR=zer os(I ength(T),I ength(f)); 

rati o_DT =zeros(I ength(T), I ength(f)); 

FR_a=zeros(l ength(T), I ength(f)); 
predl=zeros(l ength(T), I ength(f)); 
predll=zeros(l ength(T), I ength(f)); 

%% Calculate Power to acomplish temperature 

resi st=3. 675e - 10*TT+9. 93e-S; % resistivity of Platinum 

R0=resi st*l y/A; % Resistance of the filaments over temperature 

P=l 0 A 2. * RO; 

%% Calibration report 
%% input p a r a me t e r s 
% sensitivity 

suh_250Hz=122. 5; % sensitivity (high gain) at 2 5 0 H z [mV/Pa] 

sul _250Hz=0. 9732; % sensi t i vi ty (I ow gai n) at 2 5 0 Hz [mV/Pa] 

%s e n s i t i v i t y corner frequencies 
fc1u =77; % [ Hz] 

f c2u =791; % [ Hz] 

fc3u=1 8 4 3 7; % [ Hz ] 

f c 4u =60; % [ Hz] 

% phase corner frequencies 
c1u =30; % [ Hz] 

c 2 u =7 0 8; % [ Hz] 

c3u=19974; % [ Hz] 


c 4 u =60; 


% [ Hz ] 
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%% Calculations velocity e I e me n t sensitivity (High gain / Low gain) 

[su h] =s u h 250Hz./((sqrt(l+(f, A 2./fc3u A 2))).*(sqrt(l+(fclu A 2./f, A 2))).*(sqrt(l+(fc4 
u A 2. / f. A 2))). *( sqrt ( l+( f. A 2. /f c2u A 2)))); 

( su I ] =s u I 2 5 0 Hz. /((sqrt ( l+( f . A 2. / f c3u A 2) )).*( sqrt ( l+( fclu A 2. / f . A 2) )).*( sqrt ( l+( f c4 
u A 2. / f, A 2) ) ) . *( sqrt ( l+( f. A 2. /f c2u A 2) ) ) ) ; 

% 

%% Cal cul at i on of the properties of air and Castor oil 
% 

% T h e r ma I conductivity air and Castor oil 

k_a=- 6, 155e- 11, *T. A 3+6, 821 e- 8. *T. A 2 +4. 7 4 3 e- 5. *T- 7. 416e- 3; 

k_c =5. 5 5 6 e- 5. *T+0. 1 6 3 7; 


% Density at standard atmospheric pressure [kg/m A 3] 
r ho_ a =- 7. 4 9 9 e- 9. *T. A 3+l. 4 5 2 e- 5. *T. A 2 - 0. 0 1 0 2 8. *T +3. 153; 

r ho c = 1 0 0 0. / ( 1. 02714+7. 0 4 e - 4. * T T + 9. 6 6 e - 7. * TT. A 2 + 3. Oe - 9. * TT. A 3 - 4 . 91 e - 4* ( p / 100) ■ 

2. 6 3 3 e - 6 * ( p /100) . * TT- 4. 0 4 e - 9* ( p / 100) . *TT, A 2-8. 8 e- 11 * ( p/ 100) . * T T. A 3+l. 4 71 e - 

6*( p/ 100) A 2+9. 2e - 9*( p/ 100) A 2. *TT. A 2-3. 63e- 9*( p/100) A 3-l. 66e - 11*( p/ 100) A 3. *TT) ; 

% . 

% Specific Heat of air ( 4 0 0 deg C) and Castor oil ( 40 - 4 0 0 deg C) 

% Formula for the approximation of the specific heats 

c p_ a =7. 8 9 4 e - 11, *T. A 5-l. 916 e - 7. *T. A 4+l. 8 2 0 e - 4. *T. A 3-8. 417 e - 2. *T. A 2+19. 04, * T - 6 8 6. 8; 
c p_c =0. 0416. *T. A 2 - 1 5. 3 3 6. *T +2 6 8 4. 4; 

%% Cal cul ati on of partial p a r a me t e r s 

D_a=k_a. / (rho_a. *cp_a); % T h e r ma I diffusion coefficient air 

D_c=k_c. / (r ho_c. *cp_c); % T h e r ma I diffusion coefficient Castor oil 


f hc_a=D_a. *rho_a. *cp_a. / (2*pi *L*( 200e - 9*rhoc_Si Ni +100e - 9*rhoc_pt)); % corner 

f requency air f Hz] 

fhc_c=D_c.*rho_c.*cp_c./(2*pi*L*(200e-9*rhoc_SiNi+100e-9*rhoc_pt)); % corner 

frequency Castor oi I [ Hz ] 


f b a r _ L a =( 2*pi *L A 2. / D_ a) ’ *f; 
fbar_Lc=(2*pi *L A 2. / D_c) ’ *f; 
fbar_Aa=(2*pi *a_a A 2. / D_ a) 1 * f; 
f d_a =D_a. / (2* pi . *(a_a A 2)); 


% di mens i onI ess frequency at L (air) 

% di mens i onl ess frequency at L (Castor oil) 

% di mens i onI ess frequency at a (air) 

% Second corner frequency (air) 


Calculation for air with constant distance between the filaments 
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% Calculating the factor matrix befor the Bessel function 
C_a=(((v_a*a_a.*P)./((li*pi*ly).*D_a.*k_a)) , *exp(li*2*pi.*f))./fbar_Aa; % 

% . 

% N u me r a t o r of the Bessel function part 

numbes_a=l-(sqrt(li.*fbar_Aa).*besselk(l,sqrt(li.*fbar_Aa))); 

denombes_a=l - (( li . / fhc_a’ )*f). *( bessel k( 0, sqrtf li . *fbar_La) )- 
bessel k(0, sqrt( li . *fbar_Aa))); 

% . 

% The temperature difference between the filaments 

DT_a=C_a. * (numbes_a. / d e n o mb e s _ a); % DT for filaments in air 

%% Cal cul ati on of Frequency response in air 

DT a 0 =( ( ( v a* a a. * P) . / (( 2* pi * I y) . * D a.*k a) ) ) . *( I og( ( pi . *a a) . / ( 2* I y) )+0. 5 7 7 ) ; % 

DT at f=0 Hz 

n =1; 

FR_a=zeros( I engthf DT_a_0), I ength(f)); 
for n=l: I ength(T) 

FR a(n,:)=DT a 0( n) * (( 1. / s q r t ( 1 +((1. / (f he a (n). A 2 1 ))*( f. A 2 )))).*( 1./s q r t ( 1 +((1./(f d 
_a[n). A 2 1 ))*[f 7 A 2))))); 

n =n +1; 

end 

%% Calculations for Castor oil at various distances between the filaments 
FR_c=zeros( I engthf DT _ a _ 0), I ength(f)); 
n =1; 

for n =1: ns 

f bar _Ac =( 2* pi . *a_c *2. / D_c)' *f; % di mensi onl ess frequency 

at a (Castor oil) 

fd_c=D_c./(2*pi.*(a_c.^2)); % Second corner frequency 

(air) 

C_c=(((v_c*a_c.*P)./((li*pi*ly).*D_c.*k_c))'*exp(li*2*pi.*f)),/fbar_Ac; 

n u mb es_c=l - (sqrt ( li . *f bar_Ac). *bessel k( 1, sqrt{ li . *fbar_Ac))); 

denombes_c=l-( (li . / fhc_c' )*f). *( bessel k(0, sqrt( li . *fbar_Lc) )- 
bessel k(0,sqrt(li .*fbar_AcT)); 

DT_ c =C_ c. * ( numbes_c. / d e n o mb e s _ c) ; % DT for filaments in 

Castor oil 

DT_c_0=( (( v_c*a_c. *P). / ((2*pi *1 y). *D_c. *k_c))). *( I og( (pi . *a_c) . / (2*1 y)) +0. 577); 
%% Calculate the frequency response in Castor oil 
n n =1; 
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for n n = 1:I engt h(DT _ c _ 0) 


FR c ( n n, : ) =DT c Of nn) *( ( 1. / sqrt ( 1+( ( 1. / (f he c ( n n). A 2 '))*( f. A 2 )))).*( 1./s q r t ( 1 +( ( 1. / 
(fd.c(nn). A 2'))*(f. A 2))))|; 

n n =n +1; 

end 

n n =1; 

% Loop to compare performance at temperatures to performance in air 
for nn=l: I engthf DT _ a _ 0) 
r at i o _ DT ( nn, : ) =DT_ a ( 361,: ) , / DT _ c ( n n, :); 
n n =n +1; 

end 

%% Calculation of the ratio of the frequency response 
n n =1; 

% Loop to compare performance at temperatures to performance in air 

for nn=l: I ength(T) % Loop to compare performance at temperatures to performance 
i n ai r 

rat i o_FR( nn, : ) =FR_a( 361, : ). / FR_c( nn, : ); 
rati o _ F R _ c o I d (n n, : ) =FR_a( 274, : )./FR_c( nn, : ); 
n n =n n +1; 
end 
n n =1; 

for nn=l: I engthf DT _ a _ 0) 

predlfnn, :) =su_h./abs(rati o_FR( nn, :)); % predicted performance at 313 deg C 
using calibration report 

predllfnn, : )=su_h./abs(rati o_DT(nn, : )); % predicted performance at 313 deg C 
using calibration report and raw DT 

n n =n n +1; 

end 

%% Get e x p e r i me n t a I results 

d a t a =x I sreadf' C:\Users\Marni x \ Document s \ MATLAB \ Thes i s\Fi nal 
versi o n s \ Graphs! Predicti ons\test_2.xls ); 

S v =[ dat a(: , 1) d a t a (: , 2) ]; % get data on velocity sensor only 

S v ( : , 2) =10. A ( (Svf: , 2). / 20) +9); % Response in [mV/Pa] 

% . 

% Plot prediction at different distances between the wires 
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f i gu r e( 100+n) 

plot(f,predl(274,:)) % 313 degrees C 

hoi d on 

pi ot(f, predl( 161,: ), ' r ) % 2 0 0 degrees C 

s cat t er(Sv(:, 1), Sv(: , 2). /10); 

xl abel (' Frequency [ Hz]' , ' FontSi ze' , 14) 

yl abel (‘ Sensi ti vi ty [ mV /Pa] ’ , ’ FontSi ze 1 , 14) 

set(fi ndobj ( type’ , ' axes’ ), 'fontsi ze' , 12, ' FontWei ght’ , ' demi ) 

axi s([ 0 1 0 0 0 0 0.5 )) 

I egend(‘ 313 deg C',’200 deg C’,' experl me nt’, '313 deg C 100 
mi crons’ , ' Locat i on’ , ' Nort hEast' ) 

tt = strcat) {’ Predi cted sensitivity, no gain, 20% loss for }, num2str(a_c*le6),{ 
mi crons' }); 

ti tl e(tt, ' FontSi ze' , 18, ’ Font Wei ght’ , ' bol d ) 
f i gure( 200+n) 

plot(f,predll(274,:)) % 313 degrees C 

hoi d on 

pi ot(f, predll) 161, : ), ’ r’ ) % 200 degrees C 

scatter)Sv{:,1),Sv(:,2)./10); 

xl abel (’ Frequency [ Hz] ’ , ’ FontSi ze’ , 14) 

yl abel (’ Sensi ti vi ty [ mV /Pa] ’ , ’ FontSi ze’ , 14) 

set(fi ndobj ( type' , ’ axes’ ), 'fontsi ze’ , 12, ’ FontWei ght’ , ’ demi ’ ) 

axis([0 1 0 0 0 0 5 ]) 

I egend]’ 313 deg C','200 deg C’ , ’ experi ment’ , ’ 313 deg C 100 
mi crons’ , ' fontsi ze', 12, ’ Locati on’ , ’ Nort hEast' ) 

tt = strcat])' Predicted sensitivity 140dB using DT, no gain, 20% loss for }, 
n u m2 s t r (a _ c * 1 e 6), { mi c r o n s }); 

ti tl e( tt, ’ FontSi ze’ , 18, ’ FontWei ght’ , ’ bol d ) 

% . 

figure (1) % Plot ratio of DT for various distances between the 

wi r es 

piot(f,predi) 274,: )) 
axis([0 1 0 0 0 0 0.5 ]) 
hoi d al I 

xl abel (’ Frequency [ Hz] ’ , ’ FontSi ze’ , 16) 
yl abel (’ Sensi ti vi ty [ mV /Pa]’,'FontSize',16) 

66 



I e g e n d ( ' 2 0 u m' , ' 40um’ , ' 6 0 u m' , ' 8 0 u m' , ' lOOum' , ' fontsi z e' , 14, ’ Locati on’ , ' East Outsi d e' ) 

titlef'Predicted sensitivity, no gain at 313 deg 
C ’ , ' FontSi z e' , 20, 'Font We i g h t 1 , ' bold') 

figure (2) 

pi ot(f, pr ed 1 (161,:)) 
axi s([ 0 1 0 0 0 0 0.5 ]) 
hoi d al I 

xl abel (' Frequency [ Hz]' ,' FontSi ze' , 14) 
yl abel (' Sensi ti vi ty [mV/Pa]' , ’ FontSi ze' , 14) 

I egend]' 20um' , ' 40um' , ' 6 0 u m’ , ' 8 0 u m’ , 1 10 0um' ,' Locati on' , ' EastOutsi d e ’ ); 

t i 11 e(' Pr edi ct ed sensitivity, no gain at 2 0 0 d e g 
C’, ' FontSi ze', 20,’ Font We ight’,’ bold’) 

figure (3) % Plot ratio of DT for various distances between the 

wi r es 

pi ot(f, pr ed 1 1 ( 2 74,:)) 
axi s([ 0 1 0 0 0 0 0.3 ]) 
hoi d al I 

xlabel (' Frequency [ Hz]' ,' FontSi ze' , 16) 
yl abel (' Sensi ti vi ty [mV/Pa]' , ' FontSi ze' , 16) 

I egend]' 2 0 u m' , ' 4 0 u m’ , ' 6 0 u m’ , ’ 8 0 u m' , ' 100um' , ' fontsi ze' , 14, ' Locati on 1 , ' EastOutsi d e' ) 

t i 11 e(' Pr edi ct ed sensitivity using DT, no gain at 313 deg 

C,' FontSi ze', 20, 'Font We ight’, ’bold') 

figure (4) 

pi ot (f, pr ed11 (161,:)) 
axi s([ 0 1 0 0 0 0 0.3 ]) 
hoi d al I 

xl abel (' Frequency [ Hz]' ,' FontSi ze' , 14) 
yl abel (' Sensi ti vi ty [mV/Pa]' , ’ FontSi ze' , 14) 

I egend]' 2 0 u m' , ' 40um' , ' 60um’ , ' 80um' , ' lOOum' , ' Locat i on' , ' EastOutsi de' ); 

ti ti e( ‘ Predi cted sensitivity using DT, no gain at 2 0 0 d e g 

C', ' FontSi ze', 20, 'Font We ight 1 , 'bold') 

% 

figure (300 +n) 

i ma g e s c (f, T ( 1: 2 7 4 ) - 2 7 3, p r e d 1 ( 1: 2 7 4, : ) ) 
xlabel (‘ Frequency [ Hz] ’ , ' FontSi ze' , 14) 
yl abel (‘ Temperature [deg C]' , ' FontSi ze' , 14) 
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col orbar(' vert' ) 
caxis((0 .3]) 
h =coI orbar; 

yl abel (h, ’ mV/ Pa’ , ' FontSi ze’ , 14) 

tt = s t r c a t ({' Sensitivity prediction vs filament temperature for }, 
n u m2 s t r (a _ c * 1 e 6), { mi c r o n s }); 

ti tl e(tt, ’ FontSi ze' , 18, ' FontWei ght' , ' bol d' ) 

figure (400 +n) 

i ma g e s c (f, T ( 1: 2 7 4 ) ■ 2 7 3, p r e d 1 1 ( 1: 2 7 4, : )) 
xl abel (' Frequency [ Hz]' ,' FontSi ze' , 14) 
yl abel (' Temperature [deg C] ’ , ’ FontSi ze’ , 14) 
colorbar('vert') 
caxis([0 .3]) 
h =coI orbar; 

yl abel (h, ' mV/ Pa’ , ' FontSi ze’ , 14) 

tt = strcat( {' Sensi ti vi ty prediction vs filament temperature using DT for }, 
num2str(a_c*le6), { mi crons }); 

titl e(tt, ’ FontSi ze’ , 18, ’ Font We i ght' , ' bol d ) 


n =n +1; 

a_c=a_c+st ep; 
end 
n =1; 

for n =1: ns 
fi g u r e ( 100+n) 
hoi d on 

pi ot (f, pr edl( 2 7 4,; ) , ' - b' ) 

I egend(' 313 deg C','200 deg C’ , ’ experi ment' ,' 313 deg C 100 
mi crons’ , ’ Locat ion',' Nort h E a s t' ) 

n =n +1; 

end 

for n =1; ns 
f i gur e( 2 0 0 +n) 
hoi d on 


pi ot (f, pr edll( 2 7 4, :),’-- b’ ) 
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I egend( ' 313 deg C’ , ' 2 0 0 deg C" , ' experi merit' , ' 313 deg C 100 
mi crons’ , ' Locat i on’ , ' Nort hEast' ) 


n =n +1; 
end 

%% predicti ons at 313 deg Celci us 
f i gur e(5) 

pi ot (f, 10. *pr edl( 2 7 4, : )) 
hoi d on 

pi ot (10. * pr edllf 2 7 4, : ), ' r' ) 
s cat t er (Sv(:, 1), Sv(: , 2)); 

I egend(‘ predi cti on using approximation','prediction using DT ratio',' 
experi mental resul t s' , ' Locati on' , ' Nort hEast’ ) 

xl abel (' Frequency [ Hz]' , ' FontSi ze' , 16) 

yl abel (' Sensi ti vi ty [ mV /Pa]', ' FontSize' , 16) 

set(fi ndobj ( type' , ' axes' ), 'fontsi ze' , 14, ' FontWei ght' , ' demi ' ) 

axis([0 1 0 0 0 0 1 ]) 

ti tl e(' Sensi ti vi ty 14 0 d B, no gain, 20% loss','FontSize' ,20, 'Font We ight 


averaged 


, ' bol d’ ) 
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